function f=mynormcdf(x,mu,sigma)

z=(x-mu)/sigma;
f=(1/2)*erfc(-z/sqrt(2));

    %
